Orbital-order melting in rare-earth manganites: the role of super-exchange 
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We study the mechanism of orbital-order melting observed at temperature Too in the series of 
rare-earth manganites. We find that many-body super-exchange yields a transition-temperature 
Tkk that decreases with decreasing rare-earth radius, and increases with pressure, opposite to the 
experimental Too- We show that the tetragonal crystal- field splitting reduces Tkk further increasing 
the discrepancies with experiments. This proves that super-exchange effects, although very efficient, 
in the light of the experimentally observed trends, play a minor role for the melting of orbital ordering 
in rare-earth manganites. 
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The role of orbital degrees of freedom [1] in the physics 
of LaMn03, and in particular the co-operative Jahn- 
Teller transition, has been debated since long (THEj. Ab- 
initio LDA+Z7 calculations show that Coulomb repulsion 
effects are key to understanding the orbitally-ordered 
antiferro- magnetic ground state Super-exchange 
alone, however, is not sufficient Q to explain the presence 
of co-operative Jahn- Teller distortions in nanoclusters up 
to T -1150 K @, (orbitally disordered phase). Still, 
superexchange effects are rather large: Tkk, the temper- 
ature at which superexchange alone would drive the tran- 
sition, is remarkably close to Too, the temperature at 
which the co-operative Jahn- Teller distortion disappears 
in resonant X-ray and neutron scattering [9]. This fact 
could indicate that super-exchange, although insufficient 
to explain the persistence of Jahn- Teller distortions in the 
orbitally disordered phase, plays a major role in the or- 
bital order-to-disorder transition (orbital order melting) 
observed at Too- Here we resolve this issue. 

Remarkably, orbital-order melting has been observed 
[Icil-[l2| in the full series of orthorhombic rare-earth (RE) 
manganites REMnOs. These systems are perovskites 
(Fig. [1]) with electronic configuration Mn 3d 4 (t^g^g). In 
the co-operative Jahn- Teller phase (T < Too), the MnC>6 
octahedra are tilted and rotated, and exhibit a sizable 
Jahn- Teller distortion with long and short MnO bonds 
antiferro-ordered in the xy plane, and ferro-ordered along 
z (Fig. [I]). Neutron and X-ray diffraction data show that 
Too increases from 750 K to — 1500 K with decreas- 
ing ionic radius IR (La Dy) 

under increasing; 

pressure eventually orbital order melts 13|, LL4[ , while JT 
distortions still persist in nanoclusters §[. 

The strength of super-exchange is directly linked to the 
amplitude of the hopping integrals, which depend on the 
cell volume and distortions. In the REMn03 series the 
volume decreases with ionic radius. Tilting and rotation, 
however, increase, because of the increasing mismatch 
between the Mn-0 and RE-0 bond- lengths. For LaMnOs 
a volume collapse at Too has been reported [TBI- Under 
pressure, up to P=18 GPa the volume decreases by ~ 
10%, while tilting/rotation slightly decrease. A sizable 



FIG. 1: (Color online) Orbital-order in TbMnOs, as obtained 
by LDA+DMFT calculations. The pseudo-cubic axes point- 
ing along Mn-Mn bonds are shown in the left corner. 



volume reduction typically increases the Mn-0 hopping 
integrals, while tilting and rotation tends to reduce them, 
reducing super-exchange effects. The scenario is further 
complicated by the local crystal- field 0, H, [I?} , which 
can, depending on its size and symmetry, help or compete 
with super-exchange, and thus even reverse the trends. 

In this Letter we clarify the role of super-exchange in 
orbital- melting. We show that, already in the absence 
of crystal- field splitting, only in LaMnOs Tkk ^ ^bo, 
while in all other systems Tkk is 2-3 times smaller than 
Too: While Too strongly increases with decreasing IR, 
Tkk slightly decreases. Taking the tetragonal crystal- 
field splitting into account, these trends are enhanced 
even further. This proves that, although very large, in 
view of the reported experimental trends, super-exchange 
plays a minor role in the orbital-melting transition. In 
order to quantify the role of many-body super-exchange 
in determining Too, we perform ab-initio calculations 
based on the local density approximation (LP A) + dy- 
namical mean-field theory (DMFT) method [18| in the 
paramagnetic phase. The minimal model Hamiltonian 
to study super-exchange effects in manganites is the 
Hubbard model for the e g bands in the magnetic field 
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FIG. 2: Energy gain per formula unit due to orbital-ordering 
(LaMnOs). Error bars are smaller than the symbols. 



h = JSt 2g of disordered t^g spins S t2g [19j 

i i^i' mam' a' 

im im 

+ ^ ^2 ( U - J^a,cj')nimcjnim'cj' • (1) 

im(^m')aa' 



ima 



creates an electron with spin a =ff,J| in a 
Wannier orbital |ra) = |x 2 — ?/ 2 ) or 1 3^r 2 — 1) at site 
z, and n im(J = cJ m(7 c- m(T . ft (4) indicates the spin 
parallel (antiparallel) to the t<i g spins (on that site). 
In the paramagnetic state, the matrix u (uf a , = 2/3) 
accounts for the orientational disorder of the t<i g spins 
[HI; il£ m , is the LDA [io| hopping integral from or- 
bital m on site i to orbital m' on site i f ^ i, ob- 
tained ab-initio by down-folding the LDA bands and 
constructing a localized Wannier basis. The on- 
site term ajrUx + St^s yields the LDA crystal-field 
matrix. It is the sum of a Jahn- Teller {ejT^ix) and 
a tetragonal {stUz) term, where n x and Tj iZ are the 
pseudospin-1/2 operators T{ x = \^ J(T 

T iz = \ T,a,m(~ 1 ) 5rri,x2 ~ y2 4ma c irna- U and J are the di- 
rect and exchange screened on-site Coulomb interaction. 
We use the theoretical estimates J = 0.75 eV, U ~ 5 eV 
@, [HI, H3| and 2JS t2g ~ 2.7 eV [22]; we find that, in 
the high-spin regime, Tkk is not sensitive to the specific 
value of 2JSt 2g , therefore we keep h fixed in all results we 
present. We solve (pQ) within DMFT [23| using a quantum 
Monte Carlo [24] solver, working with the full self-energy 
matrix S mm / in orbital space [17j. We construct the LDA 
Wannier functions via the downfolding procedure based 
on the Nth-Order Muffin-Tin (NMTO) method Q. Ad- 
ditionally, we perform L AP W calculations [IBj , and con- 
struct maximally localized Wannier functions [26[. The 
band-structures and parameter trends obtained with the 
two methods are very similar |27j. 

To determine the super-exchange transition tempera- 
ture Tkk we use two independent approaches. In the 



first, we calculate the order parameter p as a function of 
temperature T, in the second we determine the T = 
total energy gain AE(p) (Fig. [2]) due to orbital order. 

The order parameter for orbital-ordering is the orbital 
polarization p=\ni—ri2\, where |1) and |2) are the nat- 
ural orbitals in e^-space. To determine Tkk we perform 
LDA+DMFT calculations as a function of temperature 
for all materials in the series. They differ in (i) hopping 
integrals and (ii) crystal field, due to static distortions. 
In order to separate the effects of super-exchange from 
those of the crystal field, we perform LDA+DMFT cal- 
culations of the orbital polarization as a function of tem- 
perature for the real system (iJ LDA ), for ideal structures 
with the same hopping integrals but no crystal-field split- 
ting (Fig. [3]), and for ideal structures with only tetragonal 
splitting (Fig. EJ). 

In the second approach we calculate the energy gain 
due to orbital order from the difference in total energy be- 
tween the orbitally polarized and the orbitally disordered 
states, in the absence of crystal fields (st = £jt = 0). 
We first perform LDA+DMFT calculations for decreas- 
ing temperature and calculate the total energy per for- 
mula unit and polarization p, Etot(p)- Next, we re- 
peat the same procedure, but with the constraint p = 
(Ei 5 i = E 2 ,2 and Ei j2 = 0). 

The total energy is given by [28| , 



Etot(p) = ^tot + (H) p — E t 



riLDA 



E 



DC, 



where ^ot is tne LDA total-energy, T^ DA the thermal 
average of (pQ) in the non-interacting (U = 0, J — 0) 
case, (H) p the actual thermal average of (pQ) in DMFT for 
polarization p, and Edc the double counting correction. 
Of all these terms only (H) p contributes to 

AE(p) = E TO t(p) ~ £tot(0) = (H) p - (H) p=0 . (2) 

(H) p can be split into a single-electron contribution (from 
the first three terms in (PQ)), which we calculate as sum 
on Matsubara frequencies, and a correlation contribution 
(from the last two terms in (PQ)), which we obtain from the 
double-occupancy matrix. Since — AEtot(p) ~ 10 — 50 
meV, error bars, in particular the QMC statistical er- 
ror on the double-occupancies matrix, have to be con- 
trolled to high accuracy [Isj. The total-energy gain for 
LaMnOs is shown in Fig. O We obtain similar be- 
havior for the other systems. In the zero-temperature 
limit, we extrapolate from AE(p) the super-exchange 
energy gain AE KK = E TOT (p = 1) - E TOT (p = 0). 
Remarkably, we find that the static mean-field [311 re- 
lation Tkk = |2A£?kk| Ab, which is valid for spin- 1/2 
Heisenberg-like model [l| with arbitrary coupling con- 
stants, gives transition temperatures close to those ob- 
tained from order-parameter calculations, the difference 
being a mere small shift. Our results are shown in 
Fig. El While T KK ~ T 00 in LaMn0 3 , in all other sys- 
tems Tkk is a factor 2-3 smaller than the experimen- 
tal estimate for Too- Moreover, Tkk is maximum in 
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FIG. 3: (Color online) Orbital-order transition temperature 
Tkk [13] in rare-earth manganites REMnOs versus RE 3+ ra- 
dius, with RE=Dy (triangles), Tb (squares), Nd (pentagons), 
La (circles). Full symbols: Tkk from LDA+DMFT total- 
energy. Empty symbols: Tkk from LDA+DMFT order pa- 
rameter calculations. Symbols of decreasing size: P=0 GPa, 
5.4 GPa and 9.87 GPa. Crosses: Experimental values (ambi- 
ent pressure) from Refs. [1CB12I]. 



LaMnOs, and roughly decreases with IR from RE=La 
to Tb, then increases again. Tkk also increases under 
pressure. These trends are opposite to those reported ex- 
perimentally for the orbital melting temperature. They 
can be ascribed to the increasing distortions along the 
REMnOs series, and the decrease in volume and tilt- 
ing/rotation with increasing pressure. Finally, for all 
systems super-exchange favors the occupation of the or- 
bital (signs are given for the site displayed in Fig. [4]) 
|<9) = - sin f \x 2 - y 2 ) -f cos § |3z 2 - 1), with = 90°, while 
experimentally ~ 108° in LaMnOs increasing with de- 
creasing IR to 114° in TbMn0 3 [32[. 

Due to the competition between the tetragonal crystal- 
field splitting St and super-exchange (which favor the 
occupation of different orbitals), Tkk is reduced even 
further. We find that for finite St the system is or- 
bitally ordered already at high temperature due to the 
crystal field, but the occupied orbital has = 180°. In 
Fig. [4] we show the results for St fixed at ~ 130 meV, 
sizable but smaller than for any of the considered sys- 
tems (see Fig. EJ. We find that at the reduced critical 
temperature T^, super-exchange rotates the orbital to- 
wards 90°. The change in Tkk is small for LaMnOs, but 
Tkk is reduced to 400 K for NdMnOs, and even more 
for DyMnOs and TbMnOs. Furthermore, in the zero- 
temperature limit, the smaller T^, the closer is 6 to 
180°. Thus a fixed Et ^ 130 meV enhances the trend 
found for St = 0: Tkk is larger in LaMnOs, and de- 
creases going to DyMnOs. Still, even for LaMnOs, 
is significantly larger than the experimental 108°. This 
means that a Jahn- Teller crystal- field splitting e jt is nec- 
essary to explain the experimental 0; Fig. 2] shows that 
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FIG. 4: (Color online) Rotation of the occupied state \0) as a 
function of temperature in the presence of a 130 meV tetrag- 
onal crystal field. The orbitals are shown for TbMnOs- 

such splitting has to increase for the series RE=La, Nd, 
Dy, Tb. Taking into account that tetragonal splitting 
actually increases with decreasing pressure, and substi- 
tuting La with Nd, Tb, or Dy (Fig. [5j), this trend is 
enhanced even more. For St corresponding to the real 
structures, down to 150 K we find no super-exchange 
transition for all systems but LaMnOs. These results 
can be understood qualitatively in static mean-field the- 
ory. In this approach, super-exchange yields an effective 
Jahn- Teller splitting skk = (t x )Akk, where Akk is the 
molecular field parameter; the self-consistency condition 
for orbital order is (r x ) = \ sin^tanh (^>^Je\ + £ K k/^) > 

with sin# = £kk/\/ £ t + 4k- This equation has 
a non-trivial solution (0 ^ 180°) only if Akk/2 > 
et- The critical temperature is ^kk/^kk = 
(£ T /2fe B T KK )/tanh- 1 (£ T /2fe B T KK ), with k B T° K = 
Akk/4; it decreases with increasing et-> while 9 — )• 180° 
[331 ]. For large enough et (st > Akk/2) there is no 
super-exchange driven transition at all. 

In conclusion, for the orbital-melting transition in rare- 
earth manganites REMnOs, we find that many-body 
super-exchange yields a transition temperature Tkk very 
close to Too only in LaMnOs, while in all other systems 
Tkk is less than half Too- Moreover, we find that super- 
exchange yields ~ 90° for the occupied orbital, while 
in the experimental structures ~ 108° — 114°. We also 
find that a tetragonal splitting ex reduces Tkk even fur- 
ther. €t increases substituting La with Nd, Tb or Dy 
and decreases under pressure. Finally, super-exchange 
effects become larger with increasing pressure, while ex- 
perimentally orbital order eventually melts [13, LL4| . Our 
work proves that, in the light of the experimentally ob- 
served trends, super-exchange plays a minor role in the 
orbital-melting transitions of rare-earth manganites. 

We thank I. Loa and K. Syassen for sharing unpub- 
lished data. Calculations were done on the Julich Blue 
Gene/P. We acknowledge financial support from the 



> 

CD 



Q. 

C/) 

CD 
CO 




0.5 - 



Gd Sm Pr 
DyTb Nd La 



0.5 



1.1 



1.2 



1.1 



1.2 



IR(A) 



IR(A) 



FIG. 5: (Color online) Evolution of the crystal- field (Jahn- 
Teller and tetragonal) as a function of the rare-earth ion ra- 
dius calculated for all structures reported in Refs. [ll|? E3411] . 
Filled circles of decreasing size: LaMnOs for P=0, 5.4 and 
9.87 GPa 13]. Inset: calculated occupied orbital. 
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